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Abstract 

We consider a class of Ising spin systems on a set A of sites. The sites are grouped 
into units with the property that each site belongs to either one or two units, and the 
total internal energy of the system is the sum of the energies of the individual units, 
which in turn depend only on the number of up spins in the unit. We show that under 
suitable conditions on these interactions none of the |A| Lee- Yang zeros in the complex 
z = e 2 @ h plane, where (3 is the inverse temperature and h the uniform magnetic field, 
touch the positive real axis, at least for large values of (3. In some cases one obtains, in an 
appropriately taken /3 /*■ oo limit, a gas of hard objects on a set A'; the fugacity for the 
limiting system is a rescaling of z and the Lee- Yang zeros of the new partition function also 
avoid the positive real axis. For certain forms of the energies of the individual units the 
Lee- Yang zeros of both the finite- and zero-temperature systems lie on the negative real 
axis for all (3. One zero-temperature limit of this type, for example, is a monomer-dimer 
system; our results thus generalize, to finite (3, a well-known result of Heilmann and Lieb 
that the Lee- Yang zeros of monomer-dimer systems are real and negative. 



*Dedicated to Elliott Lieb on the occasion of his eightieth birthday, in friendship and 
admiration. 



1. Introduction 

We consider a system of Ising spins on a finite set A of sites; we often think of A as a 
subset of some lattice L. Writing a = (<7i)igA) with Oi = ±1, for a spin configuration, we 
let N = N(a) = X^eA(l + a i)/2 be the total number of up spins. We will sometimes think 
of this system as a lattice gas in which r\i = (1 + 0i)/2 is the indicator of a particle at site 
i; N is then the total number of particles in the system and iV/|A|, with |A| the number 
of sites, the average density p. The average magnetization per site is m = |A| _1 ^ o~i = 
|A| _1 (2A^ — |A|) = 1p — 1. The thermodynamic properties of this system are determined 
[1] by the partition function 

Z A {z,0)= z N ^e~^-\ (1.1) 

ct:A^±1 

where U (a) is the internal energy of the configuration a, f3 is the inverse temperature, and 
z is the magnetic fugacity, that is, z = e 2/3h with h the magnetic field. Z\ is a polynomial 
in z of order |A|, with positive coefficients. 

The zeros in the complex fugacity plane of Za(z,/3), usually called Lee- Yang zeros, 
have been of interest since the original studies of Yang and Lee [2] and Lee and Yang [3] . 
For finite systems none of the |A| zeros can lie on the physically relevant positive real axis. 
But when A is a subset of some periodic lattice L and U(a) is the restriction to A of a 
translation invariant energy (with some boundary conditions), so that we may speak of the 
thermodynamic limit A / L, the zeros can in this limit approach the real axis, signaling 
(typically) the existence of a phase transition in the model [2]. The nature of the phase 
transition depends on the manner in which the zeros approach the positive z-axis as j3 or 
other parameters in U are changed. Speaking loosely, there will be a discontinuity in the 
magnetization per site, that is, a first order transition, at a value H of the magnetic field if 
the density of zeros on the real axis at z = e 2 ^ H is nonzero, and a higher order transition 
if there is a nonzero density arbitrarily close to this point [2] . 

In their second paper [3], Lee and Yang proved that for the Ising model with ferro- 
magnetic pair interactions, that is, for 

£%) = - E ( L2 ) 

{i,j}EA, i^j 

with all Jij > 0, all the zeros of Za lie on the unit circle \z\ = 1. Consequently, the only 
possible thermodynamic phase transition in this system takes place at z = 1 or h = 0. The 
Lee- Yang theorem has been extended in many ways to a variety of classical and quantum 
systems; see [4,5] for reviews. One can also prove in many cases that there is indeed a 
first order phase transition at sufficiently large /5, so that the zeros must have a nonzero 
density at z = 1 in the thermodynamic limit. 

Much less is known rigorously for general spin systems in which the zeros do not lie on 
the unit circle. This has led to numerical studies of these zeros for A a subset of a lattice L; 
in particular, the cases in which L is either 1? or the planar triangular lattice, and in which 
the internal energy is given by (1.2) with uniformly antiferromagnetic nearest-neighbor 
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interactions, that is, with J^- = JS^_j^i, J < 0, have been investigated extensively [6,7,8]. 
These systems can be proven to undergo phase transitions in the thermodynamic limit 
A / L, for large values of /3, at nonzero values of h [9,10,11]. This implies that the zeros 
of their partition functions must converge to the real axis at some point z(/3) ^ 1. More 
recently there have also been results for systems in which the zeros lie on the unit circle 
for large (3 but not for small (3 [12]. In some cases they touch the real axis, either for finite 
f3 or in the limit (3 /* oo. 

There have also been many studies of the Lee- Yang zeros of the grand canonical 
partition function for general interacting particle systems on lattices or in the continuum. 
Of particular interest to us is the case of "hard" interactions, in which for every particle 
configuration rj either U{r]) = or U{rf) = oo; put another way, some configurations are 
forbidden, while all others have no internal energy. Systems with such interactions can 
often be obtained as a suitable (3 — > oo limit of (1.2). In many interesting cases one may 
then think of the model as a system of particles (which may or may not correspond to 
the original particles) with fixed shapes, like dimers, diamonds, or hexagons, which cover 
more than one lattice site and which cannot overlap. For such systems temperature plays 
no role, so that the partition function does not depend on (3; we will write y for the 
fugacity of the new particles and Qa{v) for the corresponding partition function. It has 
been shown in particular for the case of dimers (on an arbitrary graph) that the zeros of 
Qa(u) all lie on the negative real y-axis [13]; any system with this property will of course 
not have any phase transitions in the thermodynamic limit. On the other hand, hard 
diamonds on Z 2 and hard hexagons on the triangular lattice do have a phase transition in 
the thermodynamic limit [9,14,15]. 

In this note we will first describe a new class of Ising systems for which no zeros touch 
the positive real axis, at least for large (3 (low temperature). In some of these systems all 
the zeros lie on the negative real axis, either for all values of (3 or for large (3; in others, 
the zeros are excluded from some wedge —<fi < arg z < <fi, where < < tt. We will then 
investigate the "hard" systems obtained from some of these, after suitable rescaling, in the 
limit (3 — > oo; these systems will similarly have no zeros of Q(y) encroaching on the real y- 
axis and hence no phase transition. The models obtained in this way include the monomer- 
dimer model [13] and the graph-counting models of [16,17]; our results thus generalize these 
latter results to a wider class of "hard" systems and to related low-temperature models. 

2. A class of systems with Lee- Yang zeros bounded away from the positive 
real axis 

We consider Ising spin systems decomposable into subsystems, called units, with the 
property that each site belongs to either one or two units. Examples include the (three 
dimensional) pyrochlore lattice [18,19], in which the units are tetrahedra, the (two dimen- 
sional) kagome lattice (Figure 1(a)), in which the units may be taken to be either the 
triangles or the hexagons, the checkerboard [20,12], in which the units are the alternate 
squares of the two dimensional square lattice (Figure 1(b)), and the ladder, in which ev- 
ery square is a unit (Figure 1(c)). We use the notation of Section 1 and write A a for 
the set of sites of the a th unit, with |A a | = n a , and a a for the spin configuration and 
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N a = A^ a (a a ) = N a (a) for the total number of up spins in the a th unit. Note that in 
general N(a) < J2 a N a (a) < 2iV((r), since a site i with = 1 may belong to either one 
or two units. 




(a) Kagome lattice (b) Checkerboard (c) Ladder 

Figure 1: Lattices decomposable into units with each site in two units. 

The internal energy U(a) of the system is assumed to be the sum of the internal 
energies of the units, 

U(a) = J2Ua(a a ), (2.1) 

Of 

and these are assumed to be symmetric in the spins of the unit, so that 

U a (a) = F a (N a (a)), (2.2) 

with F a a polynomial of degree at most n a . We can think of (2.2) as a mean field interaction 
among the spins in A a . The partition functions (1.1) for a single unit and for the entire 
system thus become 

rig / \ 

Z Aa (z,p) = z N ^e~^ N ^)) = £ ( n ° WW), (2.3) 

CTi=±l, ieA a 1=0 ^ ' 

Z A (z,P)= /^-^.W^). (2.4) 

<7 { =±1 

We will prove in Section 3 that under certain conditions on the function F a the Lee- 
Yang zeros of Z\(z, (3) are bounded away from the positive real z-axis at low temperature 
or, under stronger conditions, must lie on the negative real axis. The next result shows that 
such bounds on the zeros of Z\(z, (3) follow from similar bounds on the zeros of Z\ a (z, (3). 

Theorem 2.1: Suppose that the angle 4> satisfies < <p < tv/2. If each zero ( of Z\ a (z, f3) 
satisfies 

Ct^O, 7r-</><argC<7r + (/>, (2.5) 
then each zero (' of Z\(z, j3) satisfies 

CVO, tt-20< argC' < tt + 20. (2.6) 
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In particular, if each ( is real and negative, so is each 

Before proceeding to the proof of the theorem we give a simple example. 

Example 2.2: If unit a has antiferromagnetic pair interactions of equal strength between 
every pair of sites then its internal energy may be written, after adding a constant, as 

F a (N a (a ct )) = -J J2 (<Ti<T3-l) = -2\J\N a (n a -N a ), J < 0. (2.7) 

{i,j}CA a , i=£j 

We consider several cases in which each unit has energy of the form (2.7). 

(a) The one-dimensional nearest neighbor antiferromagnetic Ising model, defined on A = 
{1, . . . , n}, may be regarded as a model of this type in which the units are the pairs ctj = 
{i, i + 1}, all with the same coupling J. From (2.7) we then have Z^ ai (z, /3) = z 2 + 2az + 1 

with a = e 2 ^ J]{ > 1; this polynomial has two negative real zeros and hence, by Theorem 2.1, 
the zeros Z(z, (3) are negative real. 

(b) When the only interactions considered for the checkerboard of Figure 1(b) are pair 
interactions then one may term the system the pyrochlore checkerboard [20] , since if each 
pair of vertices in a square are connected with edges one obtains a planar representation 
of a tetrahedron. To be concrete we choose, for example, a 21 x 21 lattice with doubly 
periodic boundary conditions. It is then possible to check that with the unit energy (2.7), 
with of course n a = 4, all four zeros of Z\ a (z, (3) are on the negative real axis (this is in 
fact verified for arbitrary values of n a in Theorem 3.1 below). Theorem 2.1 then states 
that the zeros of Za will all lie on the negative real axis; since this remains true as A /'l?, 
the system will not have a phase transition at any finite temperature. In fact the pressure 
and all correlations will be analytic functions of h for all h G M. 

(c) The ladder (Figure 1(c)) illustrates the fact that two units may share several vertices 
and thus an edge. Note, however, that the form (2.2) of the total energy implies, with (2.7), 
that the coupling constant for these shared edges (vertical in Figure 1(c)) is twice that for 
the unshared (horizontal) edges. 

2.1 Proof of Theorem 2.1 

The proof of Theorem 2.1 depends on two standard results, which we quote for com- 
pleteness; see the Appendix of [21] for more details. We let A n denote the space of complex 
polynomials in z\ , . . . , z n which are separately affine in each variable, and observe that if 
P is a complex polynomial of degree at most n then there is a unique symmetric P G A n 
such that P(z, . . . , z) = P(z). A closed circular region is a closed subset K of C bounded 
by a circle or a straight line. 

Theorem 2.3 (Grace's theorem): Let P be a complex polynomial in one variable of 
degree at most n. If the n roots of P are contained in a closed circular region K and 
z x K,...,z n K, then P(zi, ...,z n )^0. 

If P is in fact of degree k with k < n then we say that n — k roots of P lie at oo and 
take K noncompact. For a proof of the result see Polya and Szego [22] V, Exercise 145. 
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Lemma 2.4 (Asano-Ruelle) [23,24]: Let K u K 2 be closed subsets of C, with K 1 ,K 2 ^ 0. 
If $ is separately affine in z\ and z 2 , and if 

$(zi, z 2 ) = A + Bz\ + Cz 2 + Dz x z 2 ^ 

whenever z\ K\ and z 2 K 2 , then 

$(z) = A + Dz^0 

whenever z £ —K\ ■ K 2 . [We have written —K\ ■ K 2 = {—uv : u G K±, v £ K 2 }]. 

The map $ i— > $ is called Asano contraction; we denote it by (z±, z 2 ) — > z. To state 
the next result we define, for e > and -tt/2 < 9 < tt/2, K e (e) = {z : Re[e l0 (z + e)] < 0}. 

Lemma 2.5: If cp is as in Theorem 2.1 and P is a complex polynomial each of whose 
zeros ( satisfies (2.5), then for any 9 with \9\ < n/2 — </> there is an e > such that 
P(z 1 ,...,z n ) t^O when z u ...,z n £ K e (e). 

Proof: Clearly there is an e > such that P(z) ^ when z ^ Kg{e), and the result follows 
from Grace's theorem. ■ 

Lemma 2.6: Suppose that <fi is as in Theorem 2.1 and that Pi(z), i = 1,...,I, is a 
polynomial of degree rii each of whose zeros ( satisfies (2.5). Suppose further that the 
polynomial Q(zi, . . . , z n ) is obtained from the product 

i 

Pi{ z i,li • • • ? z i,rii) 

i=l 

by a sequence of Asano contractions (zij,Zk,i) — > z m or relabelings Zij — > z m . Then each 
zero (' of Q(z) satisfies (2.6). 

Proof: For each 9 with |^| < n/2 — </> we obtain from Lemma 2.5 and the Asano-Ruelle 
lemma that for some e > 0, Q{z\ , . . . , z n ) ^ when each of z\ , . . . , z n lies in the complement 
of the set —Kg(e) ■ Kg(e). Thus Q(z) ^0 when z is in the complement of —Kg(e) ■ Kg(e). 
This complement is the interior of a parabola with focus at and in particular contains the 
ray making an angle 29 with the positive real axis. As 9 varies in 9 E (—tt/2 + 0, tt/2 — 0) 
this ray sweeps out the complement of the region defined by (2.6). ■ 

We can now give the proof of the main result. 

Proof of Theorem 2.1: The main statement of the theorem is an immediate consequence 
of Lemma 2.6, since Z(z\, . . . , z\\\) is obtained from T\ a Z a (z a i, . . . , z ana ) by Asano con- 
tractions and relabelings. The last statement follows by taking = 0. ■ 

In Section 3.1 below we will need the following corollary of Lemma 2.6. 

Corollary 2.7: Suppose that P(z) is a polynomial of degree n for which each zero £ is 
real and negative. If Q is obtained from P by squaring all coefficients, then each zero (' 
of Q is real and negative. 

Proof: Take = 0, / = 2 and Pi = P 2 = P in Lemma 2.6 and make all contractions 
(z ljl z 2j ) ->■ zj. m 
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3. Zeros of the partition function of a single unit 

In this section we consider a particular unit a with n a sites, energy F a (N a ), and 
partition function Z^ a , and address the question implicitly raised by Theorem 2.1: when 
are all zeros of the function Z\ a confined to a sector of the form (2.5) for some 0? In 
Section 3.1 we give a criterion which guarantees that for all (3 these zeros satisfy (2.5) with 
= 0, and in Section 3.2 several criteria implying bounds of the form (2.5) for various 
values of 0. 



3.1 A quadratic interaction energy 



Theorem 3.1: Suppose that F a (l) is quadratic with positive leading coefficient: F a {l) = 
al + bl + c with a > 0. Then for any (3 > all zeros of Z^ a (z, (3) are real and negative. 

Note that if -F a (0 = o-l 2 + bl + c then the constant c is irrelevant, the constant 
b represents a shift in the magnetic field, and the constant a may be absorbed into the 
inverse temperature; thus we may (and will) assume without loss of generality that F a {l) = 
—l(n — I). In the spin language this is an energy in which every pair of spins in the unit is 
coupled with the same antiferromagnetic interaction and there is a uniform magnetic field, 
as in Example 2.2; in the lattice gas language particles on each pair of sites interact with 
the same positive repulsive potential and there is a uniform chemical potential. 

Lemma 3.2: If b > 0, the polynomial 

P( b )(z)=^(^)(l + bl(n-l))z l 

has only real negative zeros. 
Proof: We have 



n-l 



n-2 



P (b) (z) = (z + IT + bn(n - 1) ^ ^ _ i y 

= (z + l) n + bn(n - l)z(z + l) n ~ 2 
= (z + l) n - 2 [z 2 + (2 + bn(n - l))z + 1], 
which has only real negative zeros. ■ 

Proof of Theorem 3.1: Starting from P(2~ k p) as m Lemma 3.2, we obtain by k applications 
of Corollary 2.7 that the polynomial 

j2(f)(l + mn-l)2- k fz l (3-1) 
has only real negative zeros. Letting fc-^oowe find that 

Z Aa (z^)=±(f)e^z l (3.2) 
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has only nonpositive real zeros, and we need only observe that the constant term in (3.2) 
is nonzero. ■ 



Remark 3.3: (a) If we consider the system to be comprised of a single unit, i.e., take A Q = 
A, then Theorem 3.1 implies that in the mean-field Ising model with antiferromagnetic 
interactions all Lee- Yang zeros lie on the negative real axis. 

(b) If we consider this same system but with ferromagnetic pair interactions, which is 
equivalent to taking f3 < in (3.2), then the standard Lee- Yang theory implies that all 
zeros of lie on the unit circle. 

3.2 A convex interaction energy 

When F a is as in Theorem 3.1 it is convex on the range < / < n a in the sense that 

2F a (l) <F a (l + l) + F a (l-l), l = l,...,n a -l. (3.3) 

In this section we consider a unit energy -F a (Z), not necessarily quadratic, which satisfies 

(3.3) . 

We begin by introducing some notation to describe such an F a more precisely. Let 
= ko < k\ < ■ ■ ■ < k r -i < k r = n a be indices such that strict inequality holds in (3.3) if 
and only if I = ki for some i with 1 < i < r — 1. To understand the role of these indices 
it is helpful to introduce a geometric interpretation. Let f a (x) be defined on the interval 
[0, n a ] as the linear interpolation of the nodes (/, F a (/)), / = 0, 1, . . . , n a , and let /* C ffi 2 
be the epigraph of f a : f* = {(x, y) \ x E [0, n Q ], y > f a (x)}. Then /* is a convex subset 
of 1R 2 with two vertical faces and r nonvertical faces. The vertices of /* are the nodes 
(ki,F a (ki)); all other nodes (l,F a (l)) are interior points of the (nonvertical) faces of /*. 
See Figure 2. For 1 < i < r we define H a i to be the slope of the (nonvertical) face of 
/* containing (ki-i, F a (ki-\)) and (ki,F a (ki)), and note that H ayi = F a (l) - F a (l - 1) 
whenever ki — 1 < / < ki. Finally, for h E M we define 

E a (h) = min (F a (l)-hl). (3.4) 

0<l<n 

We will be interested in Section 4 in the set S a (h) of values of / on which the minimum in 

(3.4) is realized; clearly if h is not equal to any of the then S a (h) contains a unique 
/, while for h = H a i it contains those / for which (/, F a (l)) lies in the i th nonvertical face 
of/*- 

The next result shows that at low temperature the zeros of Z\ a fall into r groups, 
where the i th group is naturally associated with the i th nonvertical face of /* and contains 
ki — ki-i points, all with magnitude of order e _/3Ha i . 

Lemma 3.4: For i = 1, . . . r let 



J — 1^1 — 1 
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Figure 2: Typical set /*, with n a = 8, r = 5, and (k , . . . , k 5 ) = (0, 3, 4, 5, 7, 8). 
Nodes (/, -F Q (/)) are indicated by dots, with heavier dots when / = ki for some i. 

and let t^i . . . U^-ki-! be the zeros of Ri. Then one may number the zeros zi, . . . z na of 
Z^ a in such a way that for ki-i < j < hi, 



converges to zero as (3 oo. Thus ki-\ of the roots converge to 0, n a — ki to infinity, and 
the remaining k{ — ki-\ to the roots of Ri [25]. ■ 

To state the main result of this section we let 5 = maxi<j< r (/cj — 5 + 1 is the 

maximum number of nodes lying on any nonvertical face of /*. 

Theorem 3.5: (a) If 5 = 1, i.e, if all the inequalities in (3.3) are strict, then all roots of 
Z\ a are real and negative for sufficiently large f3. 

(b) If 5 = 2 then all roots of Z^ a satisfy (2.5) with (p = it/ 3 for sufficiently large (3. 

(c) If 5 = 3 then there is an angle 4> na , which as indicated may be chosen to depend only 
on n a , such that 4> na < it/ 2 and such that all roots of Z^ a satisfy (2.5) with <j> = (p na for 
sufficiently large f3. 

(d) If k\ < 4 and k r -\ < n a — 4 with at least one of these an equality, and ki — < 3 
for i — 2, . . . , r — 1 so that 5 = 4, then there is an angle 4>n a < tt/2 such that all roots of 
Z\ a satisfy (2.5) with <p = 4>n a for sufficiently large (3. 

Proof: (a) When 6 = 1 each of the polynomials Ri(t) is linear, with a negative real root. 
From Lemma 3.4 the roots of Z^ a for large (3 must be widely separated in magnitude, and 
since any complex roots among these occur in complex conjugate pairs, the roots must in 
fact be real. 



lim Zje p = 



t 



i,j—ki-!- 



Proof: For / < or I > ki the coefficient of t 



in the polynomial 



Rp,i(t) = ^ EaiH ^Z Aa (te^,p) 

= t k ^Ri(t)+ Yl 




(b) It follows from Lemma 3.4 that argZj — > a,rgt i j-k i _ 1 as j3 /oo, so that it suffices to 
show that each root satisfies | argtij —n\ < tt/3. When 6 = 2 the Ri(t) are either linear, 
with roots having argument 7r, or quadratic; in the latter case the quadratic formula shows 
that the roots are complex, have negative real part, and satisfy 



Im tij 



1^6 t-ij 



which yields the desired bound. Improved bounds on the roots for specific values of ki and 
n Q may be obtained from (3.5). For example, if i = 1 then k\ = 2 and one obtains (2.5) 
with </) = 7r/4 for the smallest (in magnitude, at large f3) two zeros of Z\ a , a result closely 
related to earlier work of Ruelle [16], as we discuss in Section 5. 

(c) When 5 = 3 the polynomials Ri can be linear, quadratic, or cubic; following the analysis 
of (b) it suffices to show that in the cubic case all roots have negative real part. Up to a 
constant factor any such cubic Ri has the form 

k(k - l)(fc - 2) + k(k - l){n + 3 - k)t 

+ k(n + 3 - k)(n + 2 - k)t 2 + (n + 3 - k)(n + 2 - k)(n + 1 - k)t 3 , 

where n = n a and k = ki. For n = k = 3 this polynomial has a triple root at z = — 1; 
we may then vary n and k continuously to some desired values and ask whether roots can 
cross the imaginary axis during this procedure. We may assume that the intermediate 
values of n, k remain real and satisfy n, k > 3 and n — k > 0. Suppose (3.6) vanishes for 
t = is, s real. We cannot have s = since k(k + l)(k + 2) ^ 0. For s^Owe would have 
both 

(k - l)(k - 2) = (n - k + 2){n -k + 3)s 2 

and 

k(k-l) = (n-k+l)(n-k + 2)s 2 , 

so that (n — k + l)(k — 2) = k(n — k + 3), i.e., n + 1 = 0, in contradiction with n > 3. 

(d) We suppose that k\ = 4; the analysis when k r -\ = n a — 4 is the same. With the results 
(a)-(c) above it suffices to show that the roots of 

Mt)=p^y 0.7) 

satisfy an appropriate bound of the form (2.5). The roots of Ri(t) are all equal to —1 for 
n a = 4 and, treating n a as a continuous variable, can have the form t = is, s real, only 
for n a a root of n 2 + 9n — 4 = 0; as both the roots of this polynomial are less than 4 the 
roots of Ri for n a > 4 must all lie strictly in the left half plane. ■ 

Remark 3.6: (a) A classical result of Newton provides a converse to Theorem 3.5(a): if 
all roots of are real for some /3 then either equality holds for all / in (3.3), in which 
case F a (l) = A + Bl for some A,B, there are no interactions, and all roots of are 
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equal, or strict inequality holds for all / in (3.3). For a proof see [26], page 104. Other 
related results are contained in [27,28,29]. In particular, (a) of the theorem follows from 
the result of [27]. 

(b) By taking k ~ n a /2 one sees that there is no bound (2.5) on the roots of (3.6) which 
is uniform in n a and k and satisfies < it/2. On the other hand, one can show that such 
a uniform bound may be found both for the roots of (3.6) for fixed k and for the roots of 
(3.7). 



3.3 An example: quartic and quadratic interactions 



As an example we consider a unit with two and four spin interactions which satisfy 
spin flip symmetry. In the particle language described in Section 1 the energy is 

U a {a) = -K 2 (ViVj + ViVj) ~K A ^ {U i€ x Vi + U te x Vi) , 



l<i<j<n 



X C {l,...,n} 
1X1=4 



where rji = 1 — r]i, i.e., 

F a (l) = -K 2 



l\ fn a -l 
2 [ 2 



K 4 



We assume that K 2 and are not both zero. The convexity condition (3.3) is satisfied 
with strict inequality for all / if 

2K 2 + !)(/ - 2) + (n a - I - l)(n a - I - 2)] < 



for / 



, n — 1. This happens when 



6 na < nrg(K 2 + iK 4 ) < (p na , 
where the angles 4> n and 9 n are given by 



(3.8) 



tan# n = — — — 

(n — 2){n — 3) 



7T/2 < 6 n < 7T, 



if n is even, 



tan0 r 



(n-Vin-iy 37T/2 < n < 27T. 

if n is odd, 



(n-3) 2 ' 



See Figure 3. Under condition (3.8), Theorem 3.5(a) implies that Z\ a has its zeros on the 
negative real axis at low temperature, and hence by Theorem 2.1 so does Z\, if all units 
in the system are of this type. Note that in particular (3.8) includes the negative K 2 -axis, 
where we know from Theorem 3.1 that the zeros are real and negative at all temperatures. 
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Figure 3: In the shaded region (which extends to infinity in both 
the x and y directions) all zeros are real and negative at low temperature. 

For nonzero values of K 2 and K4 not satisfying (3.8), Remark 3.6(a) implies that the zeros 
do not lie exclusively on the real z-axis for any /3. 

If we specialize further to the case n a = 4, in which 9 n = tan -1 (—2) and <j) n = 37r/2, 
then the region in the space of interactions at which all zeros are on the negative real axis 
can be computed exactly. Since a change of temperature is equivalent to a rescaling of 
(K 2 , K 4 ) it is convenient to take (3 = 1; then this region is given by 



log 



V9 - 8a ; 





log 


"4a -3" 


>K 4 >< 






—00 





, ifK 2 >log(3/4), 
otherwise, 



where a = e K2 . See Figure 4. The computation follows that in the proof of Proposition 6 
of [12], and we omit details. 



4. Ground states and zero temperature limits 

Consider a system which is assembled from units, as described in Section 2, such that 
the energy F a for each unit is convex in the sense of (3.3). In this section we suppose 
further that each site belongs to exactly two units. 

Let us fix, for the moment, a magnetic field ho. The total energy of the system in 
spin configuration cr, including the magnetic energy, is then 

U(a) - 2h N(a) = ^[F a (iV Q (a)) - h N a (a)}. (4.1) 

a 

From (4.1) and (3.4) it follows that this energy is bounded below by Eq = Y2 a E a (ho). On 
the other hand, if we recall the definition of S^^ho) given below (3.4) we see that Eq is in 
fact the ground state energy of the system — the minimum value of (4.1) — if and only if it 
is possible to find a spin configuration g_ such that N a (a) G S a (ho) for each a. When this 
is true we say that the system is not frustrated. 
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2 " 



-3 J 

Figure 4: The case n a = 4. When /3 = 1 all zeros are real and negative 
if and only if (K2, K4) lies in the shaded region (which extends 
to infinity in both the x and y directions). 

Now we assume that our system is not frustrated and consider a zero temperature 
limit (3 f - 00 with a /3-dependent fugacity z((3) = e 2/3h< -^ such that h(/3) — > ho as (3 /• 00; 
specifically, for some A G K. we take 

h(P) = ho + ± (4.2) 

In the h-T phase plane (where T = 1//3 is the temperature) this corresponds to approaching 
(/i ,0) along a line with slope 2/A [30]. Then with y = e x the limiting partition function 
is 

Q A (y,h )= lim e^>Z A (z^)^) 

= lim r/(£) e -' 5 (^)- 2 ' , » w y- i? «) 

a_ 

where Q{ho) is the set of ground-state configurations. We are of course interested in the 
behavior of the zeros of Z A under the limiting process (4.3). If iV min (/io) = m\rY ae g^ h ^ N(a) 
and iV max (/i ) = maXo-gg^o) N(cr) then iV m j n zeros will converge to and |A| — iV max to 
00 [25], while the remaining zeros converge to the (nonzero) roots of y~ Nmin QA(y, ho). 

When for each a one has \S a (ho)\ = 1, that is, when ho 7^ H a i for any a, 2, there 
is a unique ground state configuration and Q(y,ho) is rather uninteresting. When there 
are many ground state configurations, however, they can in some cases be identified with 
configurations of "hard objects" and Q(y : H ) is then the partition function for these. 




= £ y N ^\ (4.3) 

(T€G(h ) 
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In the next example we illustrate these ideas by revisiting Example 2.2(b). In Section 5 
we describe a family of examples involving graph-counting polynomials. 

Example 4.1: We consider again Example 2.2(b): pair interactions on a 2L x 2L 
pyrochlore checkerboard with doubly periodic boundary conditions. The unit energy of 
the model is given in (2.7), with n a = 4 for all a. Since the energy for all units has the 
same form F a (l) = — 2| J\l(4 — I) we will omit the subscript a on F and similar quantities 
when no confusion can arise. The fields Hi = F(l) — F(l — 1) defined in Section 3.2 are 
H\ = -6|J|, H 2 = -2|J|, #3 = 2| J|, and H A = 6\J\. For H t < h < H t+1 (with H = -oo 
and -H5 = 00) the zero temperature limit along the line (4.2) is independent of A and the 
ground state configurations each have exactly I up spins in each unit, that is, N a = / for 
each a. If we take the T — > limit of the partition function along a line (4.2) with ho = Hi 
we obtain ground states in which both N a = / and N a = 1 — 1 are possible, with the total 
value of N controlled by the fugacity y = e x . The situation in the h-T plane is shown in 
Figure 5, with a typical line (4.2) for ho — Hi = — 6| J|. 



slope 2/A / / 
/ 

/ 



-6UI -2\J\ 



T 



h 



2\J\ 6IJ 



Figure 5: The h-T plane for the pyrochlore checkerboard with pair interactions. 



We may interpret the ground state configurations in terms of hard objects by con- 
sidering dimers on the new lattice A' obtained from A by shrinking each unit to a vertex 
and introducing an edge joining two of these points when the corresponding units share a 
site; A' is again a square lattice with certain periodic boundary conditions. An occupied 
site in A corresponds to a dimer covering the corresponding edge in A', so that a ground 
state for Hi < ho < i?z+i, with with N a = / for each a, corresponds to a dimer configu- 
ration on A' in which every vertex is covered by exactly / dimers. In a ground state with 
ho = Hi each site of A' is covered by either / or / — 1 dimers. Thus, for example, for 
ho = Hi = — 6| J\ these are monomer-dimer configurations; for ho = Hi = —2\J\ they are 
a restricted class of unbranched subgraphs [17]. In each case, Theorem 3.5(a) implies that 
all zeros of the partition function lie on the negative real axis, a result originally obtained 
for the monomer-dimer system (in much more generality) in [13]. 



5. Graph-counting polynomials 

Consider a graph G with sets V of vertices and E of edges, such that each edge 
connects a pair of distinct vertices. Note that such graphs may have several edges joining 
the same pair of vertices. We let do be the maximum vertex degree in G. A subgraph M 
of G is a graph with vertex set V and edge set contained in E; \M\ denotes the number 
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of edges of M and for any vertex v G V we let (1m (v) denote the vertex degree of v in M 
that is, number of edges of M incident on v. 

A graph-counting polynomial [16] is a polynomial 

Qc(y)= y m > 

Me{C) 

where (C) is the collection of subgraphs of G associated with some set C of nonnegative 
integers via (C) = {M \ d M (v) G C, Vv G V}. For C = {0, 1}, (C) is the class of dimer 
subgraphs of G; in this case it was shown by Heilmann and Lieb [13] that all the zeros of 
Qc are real and negative. For C = {0,1,2}, (C) is the set of unbranched subgraphs of 
G and it is shown in [17] that the zeros of Qc he in the left half plane; results for other 
choices of C are given in [16]. In this section we study Qc for C a nonempty interval 
C pq = {p,p + 1, . . . , q} of nonnegative integers. 

Given the graph G and the set C pq we introduce a statistical mechanical system of 
the type described in Section 2. In this system A = E; we may think of starting with a 
geometric realization of the graph and then putting a site of A at the center of each edge. 
For each v G V there is a unit a v which contains those sites of A which correspond to 
edges of G incident on v. Since in G each edge is incident on two vertices, this system has 
the property, assumed in Section 4, that every site belongs to exactly two units. Finally, 
we introduce a unit energy F(l), I = 0,1, ... , do(G), the same for all units, which satisfies 
(3.3) and is such that, with ki the indices defined in Section 3.2, (i) p = ki-i and q = ki for 
some i, and (ii) kj = kj-i + 1 for all j ^ i. In other words, the i th nonvertical face of the 
convex set /* associated with F (see Figure 2) contains the nodes (/, F(l)) for I = p, . . . , q, 
and all other faces contain exactly two nodes. 

Now consider the limit (4.3) with ho = Hi = (F(q) — F(p))/(q — p). The ground state 
configurations at magnetization Hi are precisely those in which N av G C pq for each v, 
and the assumption that the system is not frustrated is precisely the assumption that such 
configurations exist. Each such configuration, however, has an immediate interpretation as 
a subgraph of G: an edge e G E belongs to the subgraph if and only if the corresponding site 
in A is occupied (using the lattice gas language). With this identification the ground state 
configurations then give rise precisely to the subgraphs belonging to (C pq ), and the limiting 
partition function Q^{y,Hi) is the same as the graph-counting polynomial Qc pq (y)- 

The next result, which describes the behavior of the zeros of Qc for certain choices 
of p, q, follows immediately from Theorem 3.5. 

Theorem 5.1: Suppose that p and q are such that (C pq ) is nonempty. Then: 

(a) If q = p + 1 then all the nonzero roots of Qc pq (y) are real and negative. 

(b) If q = p + 2 then all the nonzero roots of Qc pq (y) satisfy (2.6) with = ir/3. 

(c) If q = p + 3 then there is an angle <fi, which may depend on p, q, and do(G) but not on 
the size of the graph, such that 4> < 7r/2 and such that all the nonzero roots of Qc (y) 
satisfy (2.6) with the angle <p. 

(d) If p = and q = 4 or p = do(g) — 4 and q = do(G) then there is an angle <fi < tt/2, 
which depends on^do(G) but may be chosen uniformly in the size of G, such that all the 
nonzero roots of Qc vq {y) satisfy (2.6) with the angle (p. 
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